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We present precise measurements of CP and CPT symmetry based on the full dataset of K — » tttt decays 
collected by the KTeV experiment at Fermi National Accelerator Laboratory during 1996, 1997, and 1999. 
This dataset contains about 15 million K — > tt°tt and 70 million K — » 7r+7r — decays. We measure the direct 
CP violation parameter_Re(e'/e) = (19.2 ±2.1) X 10 -4 . We find the K L -K S mass difference Am = (5265 ± 
10)xl0 6 hs" 1 and the K s lifetime r s = (89.62 ± 0.05) xl0~ 12 s. We test CPT symmetry by finding the phase 
of the indirect CP violation parameter e, <p e = (44.09 ± 1.00)°, and the difference of the relative phases between 
the CP violating and CP conserving decay amplitudes for K — » 7r + 7r~ (</>+_) and for K — » 7r°7r (</>oo)i A</> = 
(0.29 ± 0.31)°. These results are consistent with other experimental results and with CPT symmetry. 



1. Introduction 

Violation of CP symmetry occurs in the neutral 
kaon system in two different ways. The dominant ef- 
fect is the result of an asymmetry in the mixing of K° 
and K° such that Kl and Ks are not CP eigenstates. 
This effect is parameterized by e and is called indi- 
rect CP violation. The other effect, called direct CP 
violation, occurs in the K — > tttt decay process and is 
parameterized by e'. Direct CP violation affects the 
decay rates of K — > 7r + 7r~ and K — ^ 7r°7r° differently, 
so it is possible to measure the level of direct CP vio- 
lation by comparing rj^ and ^oo : 
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Measurements of tttt phase shifts |l| show that, in the 
absence of CPT violation, the phase of e' is approx- 
imately equal to that of e. Therefore, Re(e'/e) is a 
measure of direct CP violation and Im(e'/e) is a mea- 
sure of CPT violation. 

For small \e'/e\, Im(e'/e) is related to the phases of 
r/ + - and 7700 by 
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The relation of the complex parameters 77^ , 7700, e, 

and e' is illustrated in Fig. [1] using the central values 
measured by the KTeV experiment. 

Experimental results have established that Re(e'/e) 
is non-zero @, 0, EL 0. In 2003, KTeV reported 
Re(e'/e) = (20.7±2.8) x 10" 4 based on data from 1996 
and 19970. We now report the final measurement of 
Re(e'/e) from KTeV. The measurement is based on 
85 million reconstructed K — > tttt decays collected in 
1996 1997, and 1999. This full sample is two times 
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Figure 1: Diagram of CP violating kaon parameters. For 
this illustration, the parameters have the central values 
measured by KTeV and the value of e is scaled by a factor 
of 50. Although they appear distinct in this diagram, note 

that <j>^ and (f>oo are consistent with each other within 

experimental errors. 



larger than, and contains, the sample on which the 
previous results are based. We also present measure- 
ments of the kaon parameters Am and T5 , and tests 
of CPT symmetry based on measurements of Acj> and 

4>+- - 4>sw- 

For these results we have made significant improve- 
ments to the data analysis and the Monte Carlo simu- 
lation. The full dataset, including those data used in 
the previous analysis, has been reanalyzed using the 
improved reconstruction and simulation. These re- 
sults supersede the previously published results from 
KTeV. In this presentation, we will focus primarily on 
improvements to the K — > 7r°7r° analysis which have 
reduced the systematic uncertainty in Re(e'/e) relative 
to the previous KTeV result. 
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2. The KTeV Experiment 

The measurement of Re(e'/e) requires a source of 
Kl and Ks decays, and a detector to reconstruct 
the charged (tt + tt~) and neutral (7r°7r°) final states. 
The strategy of the KTeV experiment is to pro- 
duce two identical Kl beams, and then to pass one 
of the beams through a "regenerator." The beam 
that passes through the regenerator is called the re- 
generator beam, and the other beam is called the 
vacuum beam. The regenerator creates a coherent 
\Kl) + p \Ks) state, where p, the regeneration ampli- 
tude, is a physical property of the regenerator. The 
regenerator is designed such that most of the K — » tttt 
decays downstream of the regenerator are from the 
Ks component. The charged spectrometer is the pri- 
mary detector for reconstructing K — > tt + tt~ decays 
and the pure Cesium Iodide (Csl) calorimeter is used 
to reconstruct the four photons from K — > 7r°7r° de- 
cays. A Monte Carlo simulation is used to correct for 
the acceptance difference between K — > tttt decays in 
the two beams, which results from the very different 
Kl and Ks lifetimes. The measured quantities are the 
vacuum-to-regenerator "single ratios" for K — > tt + tt~ 
and K — * tt°tt° decay rates. These single ratios are 

proportional to | ^7— i / p\ 2 and / p\ 2 , and the ratio 

of these two quantities gives Re(e'/e) via Eq. [TJ 



2.1. The KTeV Detector 

The KTeV detector (Fig. (2) consists of a charged 
spectrometer to reconstruct K — > tt + tt~ decays, a 
pure Csl electromagnetic calorimeter to reconstruct 
K — > tt q tt° decays, a veto system to reduce back- 
ground, and a three-level trigger to select events. Two 
virtually identical neutral kaon beams arc incident on 
the detector; a movable active regenerator is placed in 
one of these beams to provide a coherent mixture of 
Kl and Ks- In this manner, we collect Kl — > tttt 
and Ks — > tttt decays simultaneously so that many 
systematic effects cancel in the ratios used to calculate 
Re(e'/e). 

The KTeV spectrometer consists of four square drift 
chambers and a large dipole magnet. Each drift cham- 
ber measures charged-particle positions in two planes 
each in the x and y views. The drift chamber planes 
have a hexagonal cell geometry formed by six field- 
shaping wires surrounding one sense wire, and the two 
planes in each view are offset to resolve position am- 
biguities. There are a total of 1972 sense wires in the 
four drift chambers. The magnet imparts a kick of 
412 MeV/c in the horizontal plane. The spectrome- 
ter measures the momenta of charged particles with 
an average resolution of 0.4%. The K — > tt + tt~ re- 
construction achieves a z-vertex resolution of 5-30 cm 
and a mass resolution of 1.5 MeV/c 2 . 

The Csl calorimeter measures the energies and po- 
sitions of photons from the electromagnetic decay of 
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Figure 2: The KTeV Detector 



the neutral pions in K — > 7r°7r° decays. It consists 
of 3100 pure Csl crystals viewed by photomultiplier 
tubes. The layout of the 1.9 x 1.9 m 2 calorimeter is 
shown in Fig. H There are 2232 2.5 x 2.5 cm 2 crys- 
tals in the central region, and 868 5 x 5 cm 2 crystals 
surrounding the smaller crystals. The crystals are all 
50 cm (27 radiation lengths) long. Momentum ana- 
lyzed electrons and positrons from Kl —* ft e^v de- 
cays (^3) are used to calibrate the Csl energy scale 
to 0.02%. The Csl calorimeter has an average energy 
resolution of 0.6%. The reconstructed decay vertex of 
the neutral pion is directly related to the energies and 
positions of the photons: 



Z^o — Z Cs I — 



r \2\JE\E2 



(3) 



where r\i is the transverse distance between the pho- 
tons at the Csl and E\ and E2 are the photon en- 
ergies. The K — > TT°TT a reconstruction achieves a z- 
vertex resolution of 20-30 cm and a mass resolution of 
1.5 MeV/c 2 . 

KTeV uses an extensive veto system to reject events 
coming from interactions in the regenerator and to re- 
duce background from kaon decays into non-7T7r final 
states. The veto system consists of a number of lead- 
scintillator detectors surrounding the decay region and 
the primary detectors. A three-level trigger, consist- 
ing of fast detector signals at Level 1, processing by 
custom electronics at Level 2, and a software filter at 
Level 3, is used to select events. 

2.2. Monte Carlo Simulation 

KTeV uses a Monte Carlo (MC) simulation to calcu- 
late the detector acceptance and to model background 
to the signal modes. The very different Kl and Ks 
lifetimes lead to very different z-vertex distributions 
in the vacuum and regenerator beams. Wc determine 
the detector acceptance as a function of kaon decay 
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Figure 3: Beamline view of the KTeV Csl calorimeter, 
showing the 868 larger outer crystals and the 2232 smaller 
inner crystals. Each beam hole size is 15 x 15 cm 2 and the 
two beam hole centers are separated by 0.3 m. 
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Figure 4: Data-MC comparison of fraction of energy in 
each of the 49 Csl crystals in an electron shower, (a) 
The fraction of energy in each of the 49 Csl crystals in 
an electron shower for data, (b) 2003 data/MC ratio (c) 
current data/MC ratio 



vertex and energy including the effects of geometry, 
detector response, and resolutions. The simulation of 
the detector geometry is based both on data and sur- 
vey measurements. Many aspects of the tracing and 
detector response are based on libraries created by 
GEANT [7| simulations. To help verify the accuracy of 
the MC simulation, we collect and study decay modes 
with approximately ten times higher statistics than 
the K — > irn signal samples, such as Kl — > n^e^i/ 
and K L -> 7T 7r 7r°. 

Many improvements have been made to the MC 
simulation since the previous result was published in 
2003. We have improved the simulation of electromag- 
netic showers to include the effects of incident parti- 
cle angles and to simulate the effects of wrapping and 
shims in the Csl calorimeter. The GEANT library 
used for the previous analysis was binned in energy 
and incident position; the effect of angles was approx- 
imated by shifting the incident position based on the 
angle of incidence. The shower library has now been 
expanded to include nine angles (-35 mrad to 35 mrad) 
for photons and 15 angles (-85 mrad to 85 mrad) for 
electrons. Electrons angles may be larger than pho- 
ton angles because of the momentum kick imparted by 
the analyzing magnet. Differences between the library 
angle and the desired angle are approximated by shift- 
ing the incident position. The particle energy cutoff 
applied in the GEANT shower library generation has 
been lowered from 600 keV to 50 keV for electrons; the 
photon cutoff of 50 keV is unchanged. Sixteen show- 
ers per bin have been generated. Energy deposits are 



corrected for energy lost in the 12 /jm mylar wrapping 
around the Csl crystals and the shims that are present 
between some rows of Csl crystals. 

The current Monte Carlo produces a significantly 
better simulation of electromagnetic showers in the 
Csl. Figure S] shows a data-MC comparison of the frac- 
tion of energy in each of the 49 Csl crystals in a shower 
relative to the total reconstructed shower energy for 
electrons from Kl — > TT^e^v decays. The majority 
of the energy is deposited in the central crystal since 
the Molicre radius of Csl is 3.8 cm. These particular 
plots are made for 16-32 GeV electrons with incident 
angles of 20-30 mrad, but the quality of agreement is 
similar for other energies and angles. The data-MC 
disagreement improves from up to 15% for the 2003 
MC to less than 5% for the current MC. This im- 
provement in the modeling of electromagnetic shower 
shapes leads to important reductions in the system- 
atic uncertainties associated with the reconstruction 
of photon showers from K — > 7r°7r° decays. 

We have improved the tracing of charged particles 
through the detector with more complete treatments 
of ionization energy loss, Bremsstrahlung, delta rays, 
and hadronic interactions in the drift chambers. The 
position resolution of the drift chambers was previ- 
ously treated as flat across the cell; the dependence 
of the resolution upon position within the cell is now 
included in the simulation. We have also have up- 
dated a number of parameters that go into the kaon 
propogation and decay calculations. 
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Figure 5: 7r + 7r~ invariant mass distribution for K —* 
tt + tv~ candidate events. The data distribution is shown as 
dots, the K — > n + n~("/) signal MC (MC Sig) is shown as 
a dotted histogram and the sum of signal and background 
MC is shown as a solid histogram. 

3. Data Analysis 

The K — > 7r + 7r~ analysis consists primarily of the 
reconstruction of tracks in the spectrometer. The ver- 
tices and momenta of the tracks are used to calcu- 
late kinematic quantities describing the decay. The 
K — > 7r + 7r~ invariant mass distributions for each 
beam are shown in Fig. [5] 

To reconstruct K — > 7r°7r° decays, we measure the 
energies and positions of each cluster of energy in the 
Csl. A number of corrections are then made to the 
measured particle energies based on our knowledge 
of the Csl performance and the reconstruction algo- 
rithm. The precision of the Csl energy and position 
reconstruction is crucial to the K — > 7r°7r° analysis 
and has been improved significantly since the previ- 
ous publication. 

We determine the energy deposit in each block of 
the Csl by converting the digitized information to en- 
ergies using constants for each channel that are deter- 
mined from the electron calibration. The laser cor- 
rection, which is measured using an in-situ laser and 
corrects for spill-to-spill drifts in each channel's gain, 
is applied to each block energy. We define "clusters," 
which are 7x7 arrays of small blocks or 3x3 arrays 
of large blocks, centered on a "seed block," which 
contains a local energy maximum. The cluster en- 
ergy must be corrected for a number of geometric and 
detector effects. We apply "block-level" corrections, 
which adjust the energy in each block that makes up 
the cluster, and "cluster-level" corrections, which are 
multiplicative corrections to the total cluster energy. 

The quality of the calibration and the Csl perfor- 
mance is evaluated by analzying electrons from the 
calibration sample with all corrections applied. The 
electron calibration for 1996, 1997, and 1999 is based 
on 1.5 billion total electrons. Figure [5] shows the E/p 
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Figure 6: electrons after all corrections, (a) E/p for 
1.5 x 10 9 electrons, (b) Energy resolution. The fine curve 
shows the momentum resolution function that has been 
subtracted from the E/p resolution to find the energy res- 
olution. 
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Figure 7: K — > ir ir m^o^o distributions for for data 
(dots) and signal MC (histogram) in the vacuum (left) and 
regenerator (right) beams. The sum of the background 
MC is also shown (stars). The dashed lines indicate our 
cuts. 



distribution and the energy resolution as a function of 
momentum after all corrections. The final energy res- 
olution of the calorimeter is ve/E ~ 2%/y/E © 0.4%, 
where E is in GeV. 

The position of a cluster is reconstructed by calcu- 
lating the fraction of energy in neighboring columns 
and rows of the cluster. The position algorithm uses 
a map that is based on the uniform photon illumina- 
tion across each crystal to convert these ratios to a 
position within the seed block. 

We use the cluster energies and positions along with 
the known pion mass to determine which pair of pho- 
tons is associated with which neutral pion from the 
kaon decay and to calculate the decay vertex, the 
center of energy, and the 7r°7r° invariant mass. The 
K -> 7r°7T invariant mass distributions for each beam 
are shown in Fig. [7J 

For K — > 7r°7r° decays, the z vertex is determined 
using only the positions and energies of the four pho- 
tons in the final state. Therefore, the measured z 
vertex is dependent upon the absolute energy scale 
of the Csl calorimeter. The energy scale is set using 
electrons from Kl — * 7r ± e =F ^ decays. A small resid- 
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Figure 8: Regenerator beam K — > z-vertex distribu- 

tion near the regenerator for 1999 data and Monte Carlo, 
(a) Uncorrected data, (b) Data with energy scale correc- 
tion applied. 



c 1.002 
O 

~ 1.0015 
0) 

£ 1.001 

o 

O 1.0005 
0J 

(0 1 

u 

W 0.9995 
0.999 
0.9985 
0.998 



current analysis 



2003 analysis 



60 80 100 120 140 160 

E K (GeV) 



Figure 9: Change in the final energy scale adjustment rel- 
ative to the 2003 analysis. The dashed line represents no 
data-MC mismatch. 



ual energy scale mismatch between data and Monte 
Carlo is removed by adjusting the energy scale in 
data such that the sharp edge in the z-vertex distri- 
bution at the regenerator matches between data and 
Monte Carlo as shown in Fig. [5J The final energy 
scale adjustment for 1997 data is shown as a func- 
tion of kaon energy in Fig. [5J the average size of the 
correction is ~0.04%. As a result of improvements 
to the simulation and reconstruction of clusters, the 
required energy scale adjustment is smaller and less 
dependent on kaon energy than in the previous anal- 
ysis. 

Background to the K — > tttt signal modes is sim- 
ulated using the Monte Carlo, normalized to data 
outside the signal region, and subtracted. In this 



analysis, we use decays from coherently regener- 
ated kaons only; any kaons that scattter with non- 
zero angle in the regenerator are treated as back- 
ground. Scattering background is the same for both 
K — > 7r + 7r - and K — > 7r°7r° decays so it can be iden- 
tified using the reconstructed transverse momentum 
of the decay products in charged mode; we use 
K — > 7r + 7r~ decays to tune the simulation of scat- 
tering background on which we must rely in neu- 
tral mode. Non-7T7r background is present due to 
the misidentification of high branching-ratio decay 
modes such as Kl — > w e^v, Kl — > tt^^v, and 
Kl — > 7r 7r°7r°. Background contributes less than 
0.1% of K -> tt+tt- data and about 1% of K -> 
7r°7r° data. Tables [TJ and [Til contain a summary of 
all the background fractions for each year. There 
are some variations in background levels among the 
years due to differences in trigger and veto require- 
ments. 



Source 


Vacuum Beam 
1997 1999 


Regenerator Beam 
1997 1999 


Regenerator Scattering 
Collimator Scattering 
Kl — * 7T e^v 
Kl — * n ± fi^ 1 u 


0.009% 0.008% 
0.032% 0.032% 
0.034% 0.030% 


0.073% 0.075% 
0.009% 0.008% 
0.001% 0.001% 
0.001% 0.001% 


Total Background 


0.074% 0.070% 


0.083% 0.085% 



Table I Summary of K — > it tv background levels 

After all event selection requirements are applied 
and background is subtracted, we have a total of 25 
million vacuum beam K — ► 7r + 7r~ decays and 6 mil- 
lion vacuum beam K —> 7r°7r° decays The numbers 
of events collected in each beam are summarized in 
Table [Till 

Table IIVI summarizes the systematic uncertainties 
on Re(e'/e). We describe the procedure for evaluating 
several important systematic uncertainties below. 

Acceptance: We use the Monte Carlo simulation to 
estimate the acceptance of the detector in momen- 
tum and z-vertex bins in each beam. We evaluate 
the quality of this simulation by comparing energy- 
reweighted z-vertex distributions in the vacuum beam 
between data and Monte Carlo. We fit a line to the 
data-MC ratio of the z-vertex distributions and call 
the slope of this line the acceptance "z-slope." A z- 
slope affects the value of Re(e'/e) by producing a bias 
between the regenerator and vacuum beams because 
of the very different z-vertex distributions in the two 
beams; we use the known difference of the mean z val- 
ues for the vacuum and regenerator beams along with 
the measured z-slope to evaluate the systematic error 
on Re(e'/e). 

Figure [TU] shows the measured z-slopes for the 
full Kl^tt + tt^, Kl — > 7r ± e =F zy, K L — ► 7r°7r°, and 
Kl — > 7r 7r°7r event samples. We use the tt + it~ z- 
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Source 


Vacuum Beam 
1996 1997 1999 


Regenerator Beam 
1996 1997 1999 


Regenerator Scattering 
Collimator Scattering 
Kl — ' n n n 
Photon Mispairing 
Hadronic Production 


0.288% 0.260% 0.258% 
0.102% 0.122% 0.120% 
0.444% 0.220% 0.301% 
0.007% 0.007% 0.008% 
0.002% 0.001% 


1.107% 1.092% 1.081% 
0.081% 0.093% 0.091% 
0.015% 0.006% 0.012% 
0.007% 0.008% 0.007% 
0.007% 0.007% 0.007% 


Total Background 


0.835% 0.603% 0.678% 


1.209% 1.197% 1.190% 



Table II Summary of K — > background levels. Note that photon mispairing is not subtracted from the data and is 
not included in the total background sum. 





Vacuum Beam Regenerator Beam 


K — > 7T + 7T 

K -» ttV 


25107242 43674208 
5968198 10180175 



Table III Summary of event totals after all selection cri- 
teria and background subtraction. 



Source 


Error on 

K -» 7T + 


Re(e'/e) (xlO" 4 ) 

7T- K -» 7T%° 


Trigger 


0.23 




0.20 


Csl cluster reconstruction 






0.75 


Track reconstruction 


0.22 






Selection efficiency 


0.23 




0.34 


Apertures 


0.30 




0.48 


Acceptance 


0.57 




0.48 


Background 


0.20 




1.07 


MC statistics 


0.20 




0.25 


Total 


0.81 




1.55 


Fitting 


0.31 


Total 


1.78 



Table IV Summary of systematic uncertainties in Re(e'/e) 
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Figure 10: Comparison of the vacuum beam z distribu- 
tions for data (dots) and MC (histogram). The data-to- 
MC ratios on the right are fit to a line, and the z-slopes 
(see text) are shown. All distributions are for the full data 
sample used in this analysis. 



slope to set the systematic uncertainty and measure 
the it e^v e z-slope as a crosscheck. For neutral mode, 
we use the high statistics 7r 7r 7r° mode to set the sys- 
tematic uncertainty because it has the same type of 
particles in the final state as 7r°7r° and is more sensitive 
than 7r°7r° to potential problems in the reconstruction 
due to close clusters, energy leakage at the Csl edges, 
and low photon energies. 

Energy Scale: The final energy scale adjustment 
ensures that the energy scale matches between data 
and MC at the regenerator edge, but we must check 
whether the data and MC energy scales remain 
matched for the full length of the decay volume. We 
check the energy scale at the downstream end of the 
decay region by studying the z-vertex distribution of 
7T 7r° pairs produced by hadronic interactions in the 
vacuum window in data and MC. To verify that this 



type of production has a comparable energy scale to 
K — > 7T°7r , we also study the z- vertex distribution 
of hadronic 7r°7r° pairs produced in the regenerator. 
The data-MC comparisons of reconstructed z vertex 
for these samples are shown in Fig. QT] 

To convert these shifts to an uncertainty in Re(e'/e), 
we consider a linearly varying energy scale distortion 
such that no adjustment is made at the regenerator 
edge and the z shift at the vacuum window is that 
measured by the hadronic vacuum window sample. 
The average energy scale distortion we apply is shown 
by the hatched region in Fig. [TT] We rule out energy 
scale distortions that vary non-linearly as a function 
of z vertex as they introduce data-MC discrepancies in 
other distributions. The systematic error on Re(e'/e) 
due to uncertainties in the K — * 7r°7r° energy scale is 
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Figure 11: Energy scale tests at the regenerator and vac- 
uum window. The difference between the reconstructed z 
positions for data and MC is plotted for the K — > 
regenerator , and vacuum window 7r°7r° samples. The 
solid point at the regenerator edge is the K — > 7r°7r° sam- 
ple; there is no difference between data and MC by con- 
struction. The open point at the regenerator edge is the 
average shift of the regenerator 7r°7r° samples for all three 
years. The points at the vacuum window are the shifts for 
the vacuum window samples for each year separately. The 
hatched region shows the range of data-MC shifts covered 
by the total systematic uncertainty from the energy scale. 
For reference, the data-MC shift at the vacuum window 
from the 2003 analysis is also plotted. 
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Figure 12: Comparisons of the reconstructed kaon mass 
vs z-vertex (top left), kaon energy (top right), minimum 
cluster separation (bottom left), and photon angle (bot- 
tom right) for 1999 data and MC. The values plotted are 
the difference between the reconstructed kaon mass for 
each bin and the nominal PDG kaon mass. 



0.65 x 10 4 ; this is a factor of two smaller than in the 
previous analysis. 



Energy Non-linearity: Some reconstructed quanti- 
ties in the analysis do not depend on the Csl energy 
scale, but are sensitive to energy non-linearities. To 
evaluate the effect of energy non-linearities on the re- 
construction, we study the way the reconstructed kaon 
mass varies with reconstructed kaon energy, kaon z 
vertex, minimum cluster separation, and incident pho- 
ton angle. Data-MC comparisons for these distribu- 
tions for the 1999 data sample are shown in Fig. [T2l 
To measure any bias resulting from the nonlineari- 
ties that cause the small data-MC differences seen in 
these distributions, we investigate adjustments to the 
cluster energies that improve the agreement between 
data and MC in the plot of reconstructed kaon mass 
vs kaon energy. We find that a 0.1%/ 100 GeV dis- 
tortion produces the best data-MC agreement for the 
1997 and 1999 datasets. Figure [TUl shows the improve- 
ment in data-MC agreement with this distortion ap- 
plied to 1999 data. The data-MC agreement in the 
reconstructed kaon mass as a function of kaon energy 
has been significantly improved compared to the pre- 
vious analysis in which a 0.7%/100 GeV distortion 
was required for 1997 data. 
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160 
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Figure 13: Effect of 0.1%/100 GeV distortion on Mjf vs 
Ek for 1999 data. The values plotted are the difference 
between the reconstructed kaon mass for each bin and the 
nominal PDG kaon mass. 



4. Results 

The value of Re(e'/e) and other kaon parameters 
Am, ts, <j> £ , and Im(e'/e) are determined using a fit- 
ting program. The fitting procedure is to minimize 
X 2 between background subtracted data yield and a 
prediction function. The prediction function uses the 
detector acceptance determined by the Monte Carlo 
simulation. The fits are performed in 10 GeV/c kaon 
momentum bins. To determine i?e(e'/e), we use the 
z integrated data yield; to measure the other kaon 
parameters a z-binned fit is performed. 

In the fit for Re(e'/e), the inputs are the observed 
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Figure 14: Re(e'/e) in subsets of the data sample. Each 
point is statistically independent. The dashed line indi- 
cates the value of Re(e'/e) for the full data sample. The 
97a run range includes the 1996 K — * 7r°7r° data. 



number of AT — >• tt + tt~ and K — ► 7r°7r° decays in each 
of twelve 10 GeV/c momentum bins. The kaon fluxes 
for K — > 7r + 7r~ and K — > 7r°7r° in each momentum 
bin, the regeneration parameters, and i?e(e'/e) are free 
parameters. CPT symmetry is assumed by setting the 

phases 0_| and (/>oo equal to the superweak phase. 

The final KTeV measurement of Re(e'/e) for the full 
1996, 1997, and 1999 combined dataset is: 

Re(e'/e) = [19.2 ± l.l(stat) ± 1.8(syst)} x 1(T 4 

= [19.2 ±2.1] x 1CT 4 . (4) 

The fit quality, quantified by x 2 l v — 22.9/21, is good. 

We perform several checks of our result by breaking 
the data into subsets and checking the consistency of 
the Re(e'/e) result. To check for any time dependence, 
we break the data into 11 run ranges with roughly 
equal statistics. We divide the data based on beam 
intensity, regenerator position, magnet polarity, and 
direction in which the tracks bend in the magnet. We 
check for dependence of the result on kaon momentum 
by breaking the data into twelve 10 GeV/c momentum 
bins. The Re(e'/e) results for these tests are shown in 
Figs. [T4l [T5l and [TBI We find consistent results in all 
of these subsamples. 

The combined result for 1996 and 1997 data only 
is Re(e'/e) = [20.0 ± l.7(stat)] x 10~ 4 . This result 
is consistent with the previously published KTeV re- 
sult, which is based on the same subset of data, of 
Re(e'/e) = [20.7 ± 1.5(stat)] x 10^ 4 @. 

The value of Re(e'/e) is consistent with other ex- 
perimental results. The weighted average of the 
new KTeV result with previous measurements is 
Re(e'/e) = [16.8 ± 1.4] x 10~ 4 ; see Fig. [T7J The con- 
sistency probability of these results is 13%. 

The regenerator beam decay distribution is sensi- 
tive to the kaon parameters t$, Am, (f> c , Re(e'/e), and 
im(e'/e). We measure these parameters by fitting the 
decay vertex distribution in the regenerator beam us- 
ing a single, z-binned fit. The five kaon parameters, 
ts, Am, <f> t , Re(e'/e), and 7m(e'/e) are free parame- 
ters of the z-binned fit. The fit thus provides the most 



Figure 15: Re(e'/e) consistency with beam intensity, re- 
generator position, magnet polarity, and track bend. The 
low and high intensity samples are from 1999 only and 
have average rates of ~ 1 x 10 11 protons/s and ~ 1.6 x 10 11 
protons/s, respectively. Reg-left and reg-right refer to the 
position of the regenerator beam in the detector. Mag+ 
and Mag- are the magnet polarity and in/out are the bend 
of the tracks in the magnet. In the polarity and bend sub- 
sets the K — > tt°tv° sample is common to both fits; the 
errors are estimated by taking the quadrature difference 
with the error for the full dataset. The dashed lines in- 
dicate the value of Re(e'/e) in the appropriate full data 
sample. 
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Figure 16: Re(e/e) in 10 GeV/c momentum bins. The 
dashed line indicates the value for the full data sample. 



general description of the data with no requirement of 
CPT invariance. All the systematic uncertainties are 
evaluated for the fit using a procedure identical to that 
used for the Re{e'/e) measurement, and accounting for 
correlations between the parameters. CPT invariance 
is imposed a posteriori, including the total errors of 
the parameters with their correlations, to obtain a 
precise measurement of Am and rg. 

This approach allows a self-consistent analysis of 
the data with and without CPT constraints. The re- 
sults are crosschecked with separate fits for Am and 
ts performed with CPT invariance imposed a priori. 
The two procedures agree to within 1.3cr sta t and the 
total uncertainties agree to within ~ 10%. 
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Figure 17: New world average for Re(e/e) combining re- 
sults from E73101, NA3lH, NA48[1, and KTeV. 



The results of the single z-binned fit are: 



Am | ^ 

TS |cpt 



A</> 



= [5269.9 ± 12.3] x 1(T 12 s, 

= [89.623 ± 0.047] x 10 6 h/s, 

= [43.76 ±0.64]°, 

= [44.06 ±0.68]°, 

= 4>e-4>sw = [0.40 ±0.56]°, 

= -3im(e'/e) = [0.30 ±0.35]° 
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Figure 18: A% 2 = 1 contours of total uncertaitny for (a) 
Am-Ts, (b) (4 f -Am and (c) rs-4>£- Larger ellipses cor- 
respond to the z-binned fit without CPT invariance as- 
sumption. Dashed lines correspond to <f> € = 4>sw CPT 
constraint. Smaller ellipes are obtained after applying this 
constraint. 
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The fit quality is good: X 2 /v = 425.4/(432-33). The 
correlations among these results are shown in Fig. 1 1 SI 
and Fig. [13 The value of S(f> = <f> € — <j)sw is consistent 
with zero as expected from CPT invariance in K°-K° 
mixing and the value of A</> is consistent with zero as 
expected from CPT invariance in a decay amplitude. 



5. Conclusion 

Using the full data sample of the KTeV experiment, 
we have made improved measurements of direct CP 
violation and other parameters of the neutral kaon 
system. All of these results supersede previous KTeV 
results. 

Assuming CPT invariance, we measure the direct 
CP violation parameter 

Re(e'/e) = [19.2 ± l.l(stat) ± 1.8(syst)] x 10 _4 (6) 
= [19.2 ±2.1] x 10" 4 . 

Also under the assumption of CPT invariance, we re- 
port new measurements of the Kl — Ks mass differ- 
ence and the Ks lifetime: 



Figure 19: A\ 2 = 1 contour for Re{e /e) vs Im(e'/e) as 
measured by KTeV compared to the measurement of nn 
phase shifts [lj and the CPT invariance expectation. 



Am = [5269.9 ± 12.3] x 10^ 12 s , . 

t s = [89.623 ± 0.047] x 10 6 h/s. 

To test CPT symmetry, we measure the phase dif- 
ferences 

A0 = -3Jm(e'/e) 
= [0.30 ±0.35]°, 

and 

<t> e -4>sw = [0.40 ±0.56]°. (9) 

These phase results are consistent with CPT invari- 
ance in both the decay amplitudes and K° — K° mix- 
ing. 

After decades of experimental effort, direct CP vio- 
lation in the neutral kaon system has now been mea- 
sured with an uncertainty of about 10%. Considerable 
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improvement in theoretical calculations of Re(e'/e) 
will be required to take advantage of this experimental 
precision. There is some optimism, however, that the 
next rounds of calculations using lattice gauge theory 
may approach a 10% uncertainty, making the precise 
measurements of e'/e equally precise tests of the Stan- 
dard Model. 
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